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The Gd2Sn207 pyrochlore Heisenberg antiferromagnet displays a phase transition to a four sub- 
lattice Neel ordered state at a critical temperature Tc ^ 1 K. The low-temperature state found via 
neutron scattering corresponds to that predicted by a classical model that considers nearest-neighbor 
antiferromagnetic exchange and long-range dipolar interactions. Despite the seemingly conventional 
nature of the ordered state, the specific heat Cv has been found to be described in the temperature 
range 350 mK < T < 800 mK by an anomalous power law, C„ ~ T^. A similar temperature 
dependence of Cv has also been reported for Gd2Ti207, another pyrochlore Heisenberg material. 
Such behavior is to be contrasted with the typical behavior expected for a three-dimensional 
antiferromagnet with conventional long-range order which is then generally accompanied by an 
exp(— A/T) behavior at lower temperature where anisotropy effects induce a gap A to collective 
spin excitations. Such anomalous behavior in Cv has been argued to be correlated to an unusual 
energy-dependence of the density of states which also seemingly manifests itself in low-temperature 
spin fluctuations found in muon spin relaxation experiments. In this paper, we report calculations of 
Cv that consider spin wave like excitations out of the Neel order observed in Gd2Sn2 07 via neutron 
scattering. We argue that the parametric Cv oc does not reflect the true low-energy excitations 
of Gd2Sn207. Rather, we find that the low-energy excitations of this material are antiferromagnetic 
magnons gapped by single-ion and dipolar anisotropy effects, and that the lowest temperature of 
350 mK considered in previous specific heat measurements accidentally happens to coincide with a 
crossover temperature below which magnons become thermally activated and Cv takes an exponen- 
tial form. We argue that further specific heat measurements that extend down to at least 100 mK 
are required in order to ascribe an unconventional description of magnetic excitations out of the 
ground state of Gd2Sn207 or to invalidate the standard picture of gapped excitations proposed 
herein. 



I. INTRODUCTION 

A. Persistent spin dynamics in pyrochlores 

A magnetic system with Heisenberg spins that sit on 
the vertices of a three-dimensional pyrochlore lattice of 
corner sharing tetrahedra and interact among themselves 
via nearest-neighbor antiferromagnetic exchange interac- 
tions is highly geometrically frustratedii^i^. Such a sys- 
tem is theoretically predicted to not develop conventional 
magnetic long range order at finite temperature for either 
classical^i^ or quantum spins^. As a result of this frustra- 
tion, real magnetic materials with antiferromagnetically 
coupled spins on this pyrochlore structure are highly sen- 
sitive to weak pcrturbativc interactions beyond nearest- 
neighbor exchange which dramatically affect the nature 
of the low temperature state. It is partially for this reason 
that the insulating R2M2O7 magnetic pyrochlore oxides 
have attracted such a great deal of attention in recent 
years^. Indeed, this family of materials has been found 
to display a variety of magnetic states and exotic low 
temperature behaviors that strongly depend on the spe- 
cific elements R and M considered^. 

In R2M2O7, the R site is occupied by a trivalent ion, 
such as diamagnetic Y^"*" or a magnetic rare earth ion 
(R=Gd3+, Tb3+, Dy3+, Ho3+, Er^+, Tm^+ , Yh^+) while 



the M site is occupied by a tetravalent ion which can be 
diamagnetic, such as Ti^+ or Sn*"*", or magnetic, such 
as Mo'^"*' or Mn''^+. Both the R and M sites form dis- 
tinct interpenetrating lattices of corner-shared tetrahe- 
dra and either site can be magnetic or non-magnetioii^. 
Such freedom allows for a large variety of phenomenol- 
ogy in the pyrochlore materials. This includes spin 
glasses like Y2Mo20'iii^i^iia and Tb2Mo207^, spin ices 
such as Ho2Ti207i^ii^iii and Dy2Ti207i^, conventional 
long range ordered materials like Gd2Ti20 , ^6,17, is, 19 g^j-^j 
Gd2Sn2 07>^ and even possible spin liquids as in the case 

of Tb2Ti207^>22^'S4^>S6. 

One common thread throughout these various ma- 
terials is that several experimental studies have found 
that, almost without exceptions^li^, all insulating rare- 
earth pyrochlore materials R2Ti207 and R2Sn207 dis- 
play temperature-independent spin dynamics down to 
To ^ O(IO^) mK. Indeed, residual low temperature 
dynamics has been found in pyrochlore magnetic mate- 
rials with low temperature states that range from not 
understood whatsoeve r^^'^^ to seemingly conventional 
long range orderedi^i^iii^'iSiiSS. We note in passing that 
persistent low-temperature spin dynamics has also been 
found in the Gd3Ga50i2 garnet (GGG)^i^ and in the 
SrCr8Ga40i9 kagome antiferromagnet^. We now briefly 
review the various experimentally observed behaviors of 
the R2Ti207 and R2Sn207 pyrochlore oxides. 
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Strong evidence for fluctuating spins down to ex- 
tremely low temperatures has been observed in Tb2Ti2 07 
where Ti^+ at the M site is non-magnetio2i. The rea- 
son for the failure of Tb2Ti207 to develop long-range 
magnetic order above 50 mK2^ despite a Curie- Weiss 
temperature, 0cw ^ — 14 K remains to this day largely 
unexplaine d^^'^^1^^ . Yb2Ti2 07 is perhaps just as intrigu- 
ing, with specific heat measurements revealing a sharp 
first order transition at Tc w 0.24 K^Si^, but with the 
spins not appearing static below since muon spin re- 
laxation (/iSR) and Mossbauer spectroscopy finds signif- 
icant spin dynamics down to the lowest temperature^^. 
Hence, the observed first order transition in Yb2Ti2 07 
seems rather unconventional. 

Ho2Ti207^'i^>ii and Dy2Ti207^ are frustrated 
ferromagnetsi^ and possess an extensive low-temperature 
magnetic entrop y similar to that of the common 
hexagonal Ih phase of water iceS^^. As such, the 
(Ho,Dy)2(Ti,Sn)207 materials are referred to as spin 
ice^. Theoretical and numerical studies have shown 
that the spin ice behavior originates from the long range 
nature of magnetic dipole-dipolc intcractions^i^i. Nu- 
merical Monte Carlo studies using non-local loop dynam- 
ics predict that those interactions should lead to long 
range order at low temperatures^. Yet, at variance 
with the numerical predictions, experimental studies of 
the Dy2Ti20,^ and Ho2Ti207i^i^iii have not found a 
transition to long range order down to 60 mK. In par- 
ticular, muon spin relaxation (/^SR)^ and neutron spin 
echo^ experiments find evidence for Ho^^ spin dynamics 
well below Tsi in Ho2Ti207. Interestingly, a recent neu- 
tron scattering study on Tb2Sn2 07 found a transition 
to a long-range ordered state at Tc « 0.87 K, with an 
analysis of the scattering intensity indicating that the 
observed state is a long-range spin ice stated. How- 
ever, even more recent /iSR studies find that the state 
at T < Tc in Tb2Sn2 07 remains dynamic down to the 
lowest tcmperaturo^^iii. Er2Ti207, like Tb2Sn207, was 
found via neutron scattering to display long range order 
below 1.2 K— . Yet, /iSR found persistent spin dynamics 
down to the lowest temperature—. 

Gd2Ti207 displays two consecutive transitions at 
r+ - 1 K and T" ~ 0.7 Ki^iii. Neutron scattering ex- 
perimentsi^ found that the magnetic state between T~ 
and T+ has one site out of four on a tetrahedral unit cell 
that is paramagnetic and fluctuating. At T < T^T, the 
fourth site orders, but remains much more dynamic than 
the three other sites. The microscopic mechanism giv- 
ing rise to the two experimentally observed states is still 
not understood. Here too, in Gd2Ti207, fiSK flnds con- 
siderable spin dynamics persisting down to 20 mK^i^. 
A phenomenological model for the density of states, 
g{e), has been proposed for the low-temperature state 
of Gd2Ti207. Most significantly, the proposed model 
for g{e) was shown to describe the peculiar tempera- 
ture dependence of the magnetic specific heat, Cy{T), in 
Gd2Ti207 which was found to be Cy{T) oc below T". 
Such Cy cx T^ behavior is rather unconventional. Indeed, 



in a conventional long-range ordered three-dimensional 
antiferromagnet, Cy ~ T^ down to a temperature where 
the temperature dependence turns to Cy ^ exp(— A/T) 
because of a gap A in the excitation spectrum induced 
by single ion anisotropy or anisotropic spin-spin interac- 
tions. 

In all the pyrochlore systems reviewed above, 
Tb2Ti207, Yb2Ti207, (Ho,Dy)2Ti207, Tb2Sn207, 
Er2Ti207, and Gd2Ti207, some theoretical lapses exist 
in our understanding of the equilibrium thermodynamic 
low-temperature state. Hence, it is perhaps not com- 
pletely surprising that the spin dynamics appears uncon- 
ventional in these materials with, in particular, a temper- 
ature independent /iSR spin polarization relaxation rate 
down to a baseline temperature Tq ^ 10^ mK. However, 
that tentative self-reassured standpoint is put on shaky 
ground by the fiSK and specific heat measurements on 
Gd2Sn207 that we now discuss. 



B. The case of 66280207 

It was first proposed that the aforementioned 
Gd2Ti2 07 material would be a good candidate for a clas- 
sical Heisenberg pyrochlore antiferromagnet with leading 
perturbations coming from long-range magnetic dipole- 
dipole interactions-tfi. The reason for this is that Gd'^^ is 
an S-state ion with half-filled 4f shell, hence orbital angu- 
lar moment L = 0, and spin S = 7/2. Spin anisotropy is 
therefore expected to be much smaller than for the above 
Ho, Dy and Tb based rare earth materials^. In that 
context, Gd2Sn207 should be similar to Gd2Ti207; the 
main difference being that Gd2Sn207 displays only one 
phase transition observed from a paramagnetic to a long- 
range ordered phase at Tc ^ 1 K — . Perhaps most inter- 
estingly, unlike Gd2Ti207, the experimentally observed 
long-range ordered phase in Gd2Sn207 corresponds to 
the one predicted by Palmer and Chalker for the classical 
Heisenberg pyrochlore antiferromagnet model with per- 
turbativc long-range dipolar interactions^. It is possible 
that the experimentally observed transition in Gd2Sn207 
corresponds to two very close transitions^!^ that arc not 
resolved^. 

From our perspective, Gd2Sn207 is an exemplar of 
the intriguing behavior discussed above. Yet, it offers 
itself as a crucial system to understand. The reasons 
are as follows: (i) as in Gd2Ti207, and all the ma- 
terials previously described, persistent low-temperature 
spin dynamics have been observed^ and (ii) again sim- 
ilarly to Gd2Ti207, unconventional power-law tempera- 
ture dependence of the magnetic specific heat has been 
found, specifically, Cy ~ T^. So here too, there may 
exists the possibility to relate a dynamical response and 
a bulk thermodynamic measurements to an unconven- 
tional density of states g{e) — a possible manifestation of 
the spectral down-shift that corresponds to the hallmark 
of highly frustrated systems. We see the experimental 
results on Gd2Sn207 as a crucial paradox to contend 
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FIG. 1: (Color online) The specific heat of Gd2Sn207 as a 
function of temperature, with an inset showing an enlarge- 
ment of the low temperature region from Bonville et al^ 
plotted on a logarithmic scale. The dashed line shows the re- 
sult of a relatively successful fit below 0.75 K. [Data were 
generously provided by P. Bonville]. 

with. Since the observed ordered state in Gd2Sn207 cor- 
responds to the one predicted by the model of Palmer and 
Chalker— , or a more refined model that includes Gd'^"'" 
single-ion anisotrop}»2£i^ and exchange interactions be- 
yond nearest-neighbo r^°i^°'^i , one could in principle fol- 
low the well-trodden road of solid state physics and con- 
ventional magnetism: with the Hamiltonian and con- 
sequential ground state known, identify the long wave- 
length excitations and, by second-quantizing them, calcu- 
late the low-temperature thermodynamic quantities. It 
turns out that this program was carried out in a prior 
work^ for a quantum version of a simple pyrochlore lat- 
tice model with nearest-neighbor antiferromagnetic ex- 
change plus long-range dipolar interactionsi^i^. What 
was found in Ref. |55| is that all spin wave excitations 
of the Heiscnberg pyrochlore antiferromagnet are pushed 
up in energy by the dipolar interactions and, as a result, 
all thermodynamic quantities show exponential tempera- 
ture dependence, ^ cxp(— A/T), at low temperatures^. 
The following question thus arises: 

Do the C„ - r2 results of Ref. fM] for GdjSn2 O7 
contradict the theoretical prediction of Ref. fSai l. and are 
the magnetic excitations of Gd2Sn20'j truly 
unconventional? 

This is the question that we ask, and aim to answer in 
this paper. 

To spell out the question above more specifically, we 
show in Fig. 1 the specific heat data on Gd2Sn207 re- 
produced from Ref. [23] • The behavior (dashed line 
in main panel) ranges from Tiow ~ 0.5 K to T^p ~ 0.8 K. 
The Tup is very close to the critical temperature, and one 



does not expect on general grounds (non-critical) power- 
law behaviors reflecting excitations out of the ground 
state to extend so close to the phase transition. Sec- 
ondly, the temperature TJow, when compared with the 
results of Ref. [55[ , is high compared to the temperature 
regime where we expect second-quantized (spin wave) ex- 
citations to describe this system. Finally, and this is the 
key aspect of the data that prompted the present work, 
we note that the Cv data at T < 0.5 K progressively 
droop below the dashed C„ ~ behavior. This is em- 
phasized in the inset of Fig. 1. Incidentally, we note from 
this plot that the ^ behavior does not provide a 
particularly good fit of the data for T G [350, 800] niK. 
The crux of the argument presented in this paper is that 
(i) the power law between Tiow and Tup is not a re- 
flection of the low-energy properties of Gd-based antifer- 
romagnetic pyrochlores and, most importantly, (ii) the 
behavior exhibited by below Tiow ~ 0.5 K (inset of 
Fig. 1) is a signature that the system is progressively en- 
tering a low-temperature regime characterized by expo- 
nentially activated spin excitations over a gap originating 
from both magnetic dipole-dipole interactions and single- 
ion anisotropy. We show below via calculations that ex- 
pand on the authors' previous work, Ref. [s^, that the 
specific heat data of Fig. 1 can be reasonably well de- 
scribed by such gapped magnetic excitations. In other 
words, we assert that the bulk thermodynamic proper- 
ties of Gd2Sn2 07, revealed by data like that shown in 
Fig. 1, are compatible with a conventional semi-classical 
long-range ordered phase. We suggest that specific heat 
measurements below Tiow and down to 100 mK could be 
used to confirm or disprove our proposal. The rest of the 
paper is organized as follows. We first present a model 
for exchange and dipole coupled spins on the pyrochlore 
lattice in the presence of a crystal field inducing single 
ion spin anisotropy. The Hamiltonian is decoupled via 
linear spin wave theory, and expressions for the quantum 
fluctuations and low temperature thermodynamic prop- 
erties are calculated. We then investigate the effects of 
second and third nearest neighbor magnetic exchange on 
the gap to spin wave excitations. Comparing the specific 
heat calculated in spin wave theory for Gd2Sn207 to that 
measured in the experiments of Bonville et al^, we use 
a maximum likelihood estimator to determine a set of 
further neighbor couplings which may be present in the 
material. Finally, we identify the zero temperature quan- 
tum fluctuations and present a schematic phase diagram 
of Gd2Sn207. 



II. MODEL HAMILTONIAN 

As reviewed in the previous section, the pyrochlore lat- 
tice has evinced much experimental and theoretical inter- 
est due to the large degree of geometrical frustration aris- 
ing from a structure consisting of corner sharing tetra- 
hedra. The cases of Gd2Ti2 07 and Gd2Sn207 are some- 
what special, in that the S = 7/2 Gd'^^ S-state ion should 



4 



have a relatively small intrinsic anisotropy when com- 
pared to other R2Ti207 pyrochlore oxides. It is known 
that the titanate (Gd2Ti207) has a complicated low tem- 
perature multi-fc magnetic structurei^. However, recent 
neutron scatterin g^"-^^ and electron spin resonance^ ex- 
periments performed on gadolinium stanate (Gd2Sn207) 
indicate that this material exhibits a fc = long range 
ordered state below ^ 1 K. As such, Gd2Sn207 should 
be reasonably well described by a general two-body 
spin interaction Hamiltonian which includes predominant 
isotropic magnetic exchange interactions up to at least 
third nearest neighbor and anisotropy in the form of in- 
teractions with the local crystal field as well as long range 
dipole-dipole interactionai^>^i^i^>^. 
Such a Hamiltonian can be written as 



~ Hex + Ti-dd + Ti-cf, 



(1) 



where the exchange, dipole-dipolc and crystal field terms 
are given by 

= E E J-o^^ab) SaiR') ■ S,{Rn (2a) 

i,a j,b 



"Hdd 



EE 

i,a j.b 



, [Sa{R^) ■ R^^,] [S,{R>)-R ^^ 




with the factors of 1 /2 having been included to avoid dou- 
ble counting. The various conventions used in Eqs. (|2a|) 
to (Pc)) are as follows: Sa{R^) is the spin located on one of 
N tetrahedra identified by the face centered cubic (FCC) 
Bravais lattice vector i?* and the site by one of four tetra- 
hedral sublattice vectors r^. Sa{R^) is assumed to be a 
full 0(3) operator satisfying Sa{R') ■ SaiR') = S(S + 1). 
Jab{Rab) gives the value of the isotropic Heisenberg ex- 
change interaction between two spins separated by iJj'j = 



R^ 



rh ~ R^ — Va, with a negative sign corresponding 



to antiferromagnetic interactions. In this study, we fo- 
cus on the Gd2Sn207 material, and thus consider a fixed 
value of Ji = 3ecw/[zS{S + 1)] = -0.273 K where the 
Curie- Weiss temperature is Qcw = —8.6 K and z = 6 
is the coordination number on the pyrochlore lattice^l. 
We treat the exchange interactions J2 and J31 beyond 
nearest neighbors as parameters to be adjusted below 
to produce agreement with experimental (specific heat) 
measurements on Gd2Sn207. Following the approach 
of Wills et al^ we treat the two possible third near- 
est neighbor (NN) exchange paths J31 and J32, known to 
be present in the pyrochlore a^°'^^ , separately (see Fig.|3]). 
With the expectation that J32 ^ J31, we henceforth set 
J32 = 0. The strength of the dipole interaction is given 
by -Ddd = /^o(ffMB)^/47''- At nearest neighbor distance, 
i?nn = aV2/4 = 3.695 A , where a = 10.45 A is the size 




FIG. 2: (Color online) A schematic portion of the pyrochlore 
lattice, detailing the paths which correspond to second nearest 
neighbor and two types of third nearest neighbor exchange 
interactions. 



of the cubic unit cell, Ddd/ Rnn' is approximately 15% of 
the exchange energy Ji . The crystal field Hamiltonian is 
written as an expansion of Stevens operators, O™, that 
transform like the real tesseral harmonica^. The number 
of terms in the expansion is strongly constrained by sym- 
metry and, from recent electron spin resonance (ESR) 
measurements^, the values of B2 and have been esti- 
mated at (47±1) mK and (0.05±0.02) mK, respectively. 
Here, we only consider the dominant lowest order term 
in the expansion of Ticf, which contributes energet- 
ically on equal footing with the dipole interactions, and 
leave the inclusion of higher order corrections to a future 
study. Writing the Steven's operators in terms of angu- 
lar momentum operators^, the crystal field part of the 
Hamiltonian 7i, Ticf , is: 



Hcf = ~^NBlS{S + 1) + 3B° V [Sa{R') 



(3) 



where the four unit vectors Za describe the local (111) 
direction for each site on a tetrahedron. The conven- 
tions and definitions used in this study for all vectors 
and lengths on the pyrochlore lattice are given in Table 1 
of a previous work by one of the authors^. 

We are interested in the effects of the low energy excita- 
tions (spin waves) on the thermodynamic properties of a 
real material described by Eq. ([T]). At zero temperature, 
we assume that the system is in one of the six discrete 
Palmer- Chalker (PC) A; = ground states^ depicted in 
Fig. [3] We have confirmed by direct numerical simula- 
tions that the classical zero temperature ground state in 
a model with nearest-nearest neighbor antiferromagnetic 
exchange and long-range dipolar interactions at the level 
of 10 - 20% of the exchange is the PC ground state. 
See also Refs. dOiHiH. The use of the PC state for 
Gd2Sn2 07 is supported by recent powder neutron scat- 
tering experiments^ where the magnetic diffraction pat- 
tern was compared to the expected result from multiple 
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FIG. 3: (Color online) Six degenerate Palmer- Chalker (PC) 
ground states (reverse spins for other three) for spins on a 
single tetrahedron. The ground states are characterized by 
each spin being parallel to an edge of the tetrahedron that 
does not intersect its vertex and all are tangent to the sphere 
which circumscribes it. The net magnetic moment on each 
tetrahedron is identically zero. 



candidate ground states. Our approach below will be 
to analyze the stability of this ground state and its ac- 
companying excitations by investigating the role of quan- 
tum fluctuations in reducing the fully polarized classical 
spin value of S'ci = (0,0, S). This can be accomplished 
by changing the axis of quantization from the global z- 
direction (an arbitrary choice) to a local axis described 
by a triad of unit vectors This triad is defined such 

that the locally quantized spin, denoted by a tilde, is re- 
lated to the spin operator in the Cartesian lab frame via 
a rotation 



(4) 



where u runs over the Cartesian indices {x, y, z} and 
points in the direction of each classical spin. Using 
Eq. ([4]) we can rewrite Eq. ([1]) in a much more compact 
form 

^ = -^EEE^"(i20JaT(J2^b)^6(J?^), (5) 

a,b i.j u.v 

where we have neglected a constant term and 



D. 



dd 



K-nl-i{<-Rl){nl-I^l) 



&Bl{n1 ■ Za)[nl ■ Zb)5a,b5u,v 



(6) 



Having manipulated our Hamiltonian into a more man- 
ageable form, we employ in the next section the methods 
of linear spin wave theory to diagonalize Eq. ([5]) . 



III. LINEAR SPIN WAVE THEORY 

In a previous study^, which we henceforth refer to 
as DG, we presented the diagonalization of Ticx + '^dd 



on the pyrochlore lattice via a Holstein-PrimakofE^ spin 
wave expansion to order 1/S, through the introduction of 
bosonic spin deviation (magnon) creation (annihilation) 
operators c\ (c^). The Ewald summation technique^ 
was used to calculate the Fourier transform of the infinite 
range dipole-dipole interaction matrix. The calculations 
of DG can be straightforwardly generalized to include 
the effects of the crystal field by shifting the diagonal 
spin interaction matrix elements {l^aa{k) and ^aa{k) 
of Eq. (16) in DG) by a term proportional to B2 (see 
Eq. ([6|)). The result, after the usual Bogoliubov diag- 
onalization procedure, is a Bose gas of non-interacting 
spin waves described by 



7^ = W(0)+}_^}^£(fe) 

k a 

where the summation is over all wavevectors in the first 
Brillouin zone (BZ) of the FCC lattice. The disper- 
sion relations for the spin wave modes, £a(fe), are cal- 
culated from the spectrum of the Bogoliubov transfor- 
mation. Physically, they arc identical to £a(fc) = huja{k) 
where uja{k) are the classical excitation frequencies ob- 
tained by linearizing the classical equations of motion for 
interacting magnetic dipoles or rotors^. 

In a real magnet, spin wave fluctuations with disper- 
sion ea{k) raise the classical ground state energy and 
reduce the staggered magnetic moment per spin from its 
classical value of S. From Eq. ([7]), the contribution to 
the ground state energy is given by 



(7) 



(8) 



The full spectrum of the Holstcin-Primakoff transforma- 
tion can be used to calculate the reduction in the stag- 
gered magnetization 



AS = 



\ k / 



(9) 



where Q is the 8x8 hyperbolically normalized matrix 
of eigenvectors such that Tr (Q^^HQ) = J2a^aik) and H 
is the 8x8 block matrix Hamiltonian [sec DG Eqs. (16) 
and (19)]. N is the number of tetrahedra on a pyrochlore 
lattice with periodic boundary conditions. 

At low temperatures (fc^T < ea{k)) expressions for the 
specific heat at constant volume, Cv, and staggered mag- 
netization, m = S — AS, can be derived from the classical 
partition function Z = Tr [exp(— /37i)] corresponding to 
Eq. dZl). Using Eqs. d?]) to (©, we find (sec DG) 



= P^J2Y.^^-('')''Biea{k))rcMf3£aik)],m 

k a 

m = S + i 

^^EE[Q^Q]-[i + "'^(^'^W)]' (11) 

k a 

where (3 is the inverse temperature and nB{£a{k)) = 
2/(g/3eo(fc) _ 2) is the Bose distribution function. 
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FIG. 4: (Color online) The first Brillouin zone of the py- 
rochlore lattice, showing the five high symmetry points F = 
(0,0,0), X = 27r/a(l,0,0), W = 27r/a(l, 1/2, 0), L = 
27r/a(l/2, 1/2, 1/2) and K = 27r/a(3/4, 3/4, 0) and the path 
in fc-space along which spectra are plotted in this study. 



IV. RESULTS 

Having determined the expressions for the zero tem- 
perature quantum fluctuations, as well as finite temper- 
ature thermodynamic relations of a dipolar Heisenberg 
model with crystal field interactions on the pyrochlore 
lattice, we may now study the quantitative effects of per- 
turbative J2 and J31 for fixed Ji and The limit of 
stability of the proposed PC classical ground states on 
the pyrochlore lattice can be investigated by searching 
for soft modes {Saik) 0) at some wavevector k. 



A. Second and third NN exchange 

If we ignore all interactions except nearest neighbor 
isotropic antiferromagnetic exchange, it has been known 
for some time^ that both classical^ and quantun*^ spins 
on the pyrochlore lattice fail to develop conventional long 
range magnetic order at nonzero temperature. There 
exists two unbroken symmetries in the ground state 
manifoldiii^ corresponding to lattice translation and time 
reversal. The excitation spectrum consists of two sets of 
two degenerate modes, the first being soft over the entire 
Brillouin zone, (see Fig. |4|) and the second are acoustic 
and linear along T ^ X leading to divergent quantum 
fluctuations. The inclusion of both dipole-dipole, and 
crystal field interactions break the rotational symmetry 
and lift all excitations to finite frequency^. 

Fig. O shows the dispersion spectrum of the four spin 
wave modes for three values of J2 and J31 along the 
high symmetry path described in Fig. U] We do indeed 
observe only optical modes, with degeneracy preserved 
along X ^ W lor J2 ~ J31 = 0. The qualitative be- 
havior of ea{k) with varying J2 and J31 also seems to 
confirm the naive expectation that even perturbatively 



4.0 - 



■3.0 



2.0 



1.0 



0.0 





















/ / 


//\ 
■■ / ' 
■7 ' 

■7 ' 


> 

\ 

N 


/\ 

\. '■^ 
/\ 

\"--. 


^\ \. 

.V 

..- \ 


/ / 

/-/ 

1 \ 


/' 
// 
li 
7' 

■1 








\\. 
\\, 
xV. 










\^ 


— ■hi 


Ji|= 0.00, J31/IJ1 


= 0.00 




--■•Vl^il = 0.03, J31/I./1 


= -0.03 




./2/I./1I = -0.03, J31/IJ1I = 0.03 





X 



w 



L 



K 



FIG. 5: (Color online) The spin wave excitation spectrum in 
kelvin for ( J2/I Ji |, J31/! Ji |) equal to (0.0,0.0), (0.03,-0.03) 
and (—0.03,0.03) plotted along a high-symmetry path in the 
first Brillouin zone of the FCC lattice. For the three parame- 
ter sets shown there exists a finite gap to spin wave excitations 
throughout the zone. 



small third nearest antiferromagnetic neighbor exchange 
should reduce the stability of the PC states^^i^i^. It 
also appears that the minimum excitation gap does not 
always occur at fe = indicating the underlying presence 
of a finite wavevector instability as J2 and J31 are tuned 
away from J2 = 1/31 = 0. 

The spin wave energy gap can be analyzed more quan- 
titatively by defining 



A(fc) 
A 



min [eQ(fc)] 

a 

min [A(A;)] . 

k 



(12 a) 
(12b) 



The value of A(fc) can be investigated as a function of J2 
and J31 at each of the high symmetry points described 
above. As we vary J2 and J31 through some critical 
values, instabilities first appear at these wavevectors of 
high symmetry. The resulting gap values are shown in 
Fig. El Although A(r) > for all values of J2 and J31 
studied here, the region of stability of the PC states is 
defined by the observed appearance of soft modes, at 
k = K = 27r/a(3/4, 3/4, 0) for ferromagnetic (positive) 
J2 and antiferromagnetic (negative) J31. Performing a 
search for the minimum value of the gap over the entire 
Brillouin zone (A) at each value of J2 and J31 confirms 
that the instability first appears at the if-point. The 
effect of perturbative second and third NN exchange in- 
teractions on the global minimum energy gap (Eq. (jl2bp ) 
along with the corresponding magnitude of spin fluctua- 
tions AS/S (Eq. ^) is shown in Fig. [71 Here we observe 
that upon reaching a pair of critical values for J2/IJ1I 
and J31/IJ1I, the excitation gap is suppressed to zero 
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FIG. 6: (Color online) The magnitude of the lowest spin wave 
excitation energy A(A;) in kelvin at various high symmetry 
points in the first BZ plotted in the J2 — J31 plane. All panels 
are plotted with the same color scale making the slight vari- 
ations in the gap at F = (0, 0, 0) difficult to discern. A soft 
mode instability occurs only for antiferromagnetic third NN 
couphng J31 
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FIG. 7: (Color online) The spin wave excitation gap in kelvin 
(A, top panel) found by determining the minimum energy 
over 19^ discrete points in the first BZ, and the reduction 
in sublattice magnetization (AS/S, bottom panel) for various 
values of J31/IJ1I plotted against J2/IJ1I. The jump in AS/S 
occurs once the limit of stability of the fc = Palmer-Chalker 
ground state is reached. Both panels share a common legend. 



(top panel), and divergent spin fluctuations ensue (bot- 
tom panel). The values of J2 and J31 corresponding to 
A ^ can be identified, and are best described by the 
linear relationship J31 = 0.750J2 — 0.077| Ji|. This hne 
defines the phase boundary between a sector of stability 
for the fc = PC ground states, and a region charac- 
terized by instabilities at finite wavevector. In addition, 
this line corresponds to the white regions in Fig. [5] where 
A — > 0, and thus defines the limit of applicability of the 
spin wave calculation around the PC ground state de- 
scribed in Section Hill 

Plotting A(A;) along T^X^W^L^K^T with 
J31 pinned to this phase boundary leads to the spectrum 
shown in Fig. [S] It is apparent from this result, that 
once the value of the third NN exchange constant has 
been set at a suitably antiferromagnetic value, altering 
the second NN exchange constant, has a relatively limited 
effect on the gap and on the consequential proliferation 
of quantum fluctuations about the classical ground state. 

To summarize, the effects of perturbative second and 
third NN exchange interactions on the appearance of soft 



modes and their accompanying quantum fluctuations in 
a model of a dipolar coupled antiferromagnetic Hcisen- 
berg pyrochlore with single-ion anisotropy is globally il- 
lustrated in Figs.[5][8l Such a model should well charac- 
terize the low temperature behavior of Gd2Sn2 07, and we 
next apply these tools with the goal of searching for the 
unconventional spin excitations believed to be present in 
this material on the basis of the unconventional Cy oc 
specific heat in the temperature range [350, 800] mK. 



B. The case of Gd2Sn2 07 

As described in the Introduction, recent studies^li^ 
of Gd2Sn207 have reported, on the basis of yuSR mea- 
surements, evidence for Gd'^''' spin dynamics well below 
0.9 K, as well as suggesting that the low temperature 
specific heat is accurately described by an anomalous 
power law. This is in stark contrast with the expected 
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FIG. 8: (Color online) The lowest excitation energy A{k) 
in kelvin plotted along the high symmetry path discussed in 
the text with (J2, J31) fixed to reside on the phase bound- 
ary parameterized by J31 = 0.750J2 — 0.077. The loca- 
tion of the zero modes in fc-space indicate that the instabil- 
ity out of the PC states occurs via excitations described by 
fc = 27r/a(3/4,3/4, 0). 
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FIG. 9: (Color online) The low temperature specific heat from 
Ref. [231 plotted against inverse temperature on a logarithmic 
scale. The lines are the calculated value of the specific heat 
using the ten best fit values for (J2/IJ1I, J31/IJ1I) (as seen in 
the inset). 



behavior for a three dimensional antiferromagnet, with 
possible exponential suppression at a temperature below 
a characteristic excitation gap. 

On the other hand, the long-range ordered state found 
by neutron scattering in Gd2Sn207 is that predicted 
by the simple model of Eq^H in Section III and dis- 
cussed in Refs. |lal4^i5Gl5llf52| and, consequently, the 
low-temperature behavior of this material should be well 
described by linear spin wave theory. Thus, in an attempt 
to resolve the paradox offered by the C„ ~ behavior, 
we have calculated the low temperature specific heat via 
Eq (jlOp within the J2 — J31 plane, and have performed 
a search for the parameters which best reproduce the re- 
ported low temperature specific heali^l. This was accom- 
plished by performing least squares linear fits of log Cy 
vs 1/T for T < 0.5 K between the experimental data and 
the spin wave specific heat for approximately 500 values 
of J2 and J31 . A characteristic subset of the large number 
of performed fits are displayed in Fig. [9l The values of J2 
and J31 which provided the best fit to the experimental 
data can be quantified by defining a maximum likelihood 
estimator which is shown in Fig. (TUl It is important to 
note that the fits of the specific heat, Cy, discussed here, 
were done with an absolute dimensionfuU scale, and thus 
no vertical adjustment of the experimental data was al- 
lowed. The comparisons allow only for adjustments of J2 
and J31 which arc therefore fine-tuning effects. As such, 
it appears that a model which possesses solely nearest- 
neighbor exchange, long-range dipolar interactions and 
single-ion anisotropy already leads to a reasonable semi- 
quantitative description of the Cy data below 500 niK. 
This indicates that the temperature T ^ 500 mK corrc- 
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FIG. 10: (Color online) The value of the maximum likelihood 
estimator (x^) in the J2 — J31 plane for a fit of the logarithm of 
the specific heat at low temperatures {T < 0.5 K) calculated 
using linear spin wave theory and compared with the exper- 
imental data of Ref. [53]. The parameters ( J2/I Ji |, J31/I Ji |) 
with the smallest x^j giving the most effective fit, fall along 
the indicated line J31 = O.76OJ2 - 0.0141 Ji 1. 



sponds to the upper temperature below which magnetic 
excitations become thermally activated. 

The minimum of in the J2 — J31 plane falls along the 
straight line J31 = O.76OJ2 - 0.014| Ji|. This line of best 
fit also falls in a region of large stability (highly gapped 
spin wave excitations) for the classical PC ground states 
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FIG. 11: (Color online) The low temperature specific heat 
(data from Ref. [S^l) versus temperature on a logarithmic 
scale. The calculated specific heat is plotted for the param- 
eter set J2 = 0.02, J31 = 0.0 which had the smallest x^i 
although any values close to the line mentioned in the text 
give a qualitatively very similar result. The agreement is quite 
good at low temperatures, where spin wave theory should be 
applicable. The low temperature suppression of C„ is char- 
acteristic of gapped excitations, in contrast with previously 
reported power law behavio j^^'^^ . The inset displays the tem- 
perature dependence of the Gd'^"'" moment calculated using a 
Bose gas of excitations along with the value measured from 
^^^Gd Mossbauer measurements^. 

with respect to quantum fluctuations (see Fig. [5]). 

The poorness of fit for simultaneously strong ferromag- 
netic second NN and antiferromagnetic third NN inter- 
actions or vice versa (top left, or lower right of Fig. [TO]), 
seems to indicate that is quite unUkcly that Gd2Sn207 
resides in these portions of the phase diagram. The pa- 
rameters J2 = 0.02 and J31 = 0.0 provide the hest empir- 
ical fit to the experimental specific heat data, although 
qualitatively similar fits are seen for all parameters which 
satisfy J31 = O.76OJ2 — 0.014| Ji|. Setting the parameters 
to these particular values, we display the spin wave spe- 
cific heat, as well as the temperature dependence of the 
order parameter m in Fig. 1111 Again, as shown in Fig. [9l 
it is clear that at low temperatures, (T < 370 mK). the 
experimental specific heat data systematically falls be- 
low the gapped spin wave results. This behavior is ten- 
tatively consistent with the fact that spin wave theory 
produces a smaller value for the magnetization m than 
what is measured from Mossbauer experiments. How- 
ever, we note that due to the intrinsic short dynamical 
time scale probed by Mossbauer measurements, the ex- 
perimental data in the inset of Fig. [TT] may not reflect 
the true value of the infinite-time order parameter. 

It is perhaps worthwhile to make a few comments on 
the physical meaning of the above fits. Firstly, we note 
that because of the weakly dispersive nature of the two 



lowest lying gapped magnon excitations (see Fig. [5]), the 
temperature dependence of thermodynamics quantities 
in the pyrochlore Heisenberg antiferromagnet plus dipo- 
lar interactions do not display the typical ^ be- 
havior as the temperature reaches approximately 0.5 K 
and exits its characteristic low-tempcrature exponential 
behavior. In fact, such an observation was already made 
in Ref. [ssj independently of any attempt to describe 
for Gd2Sn207. Secondly, the calculations presented here 
constitute a standard procedure for a system with con- 
ventional long range magnetic order. In this context, it 
is therefore interesting to note that, contrary to the re- 
ported behavior, the experimental specific heat data 
are not only relatively well fit using the exponential spin 
wave form of Eq. (fTO|) , but it appears to fall off even faster 
than the exponentials considered at low temperatures. 
This would lend credence to the view that analyzing ex- 
perimental data on a log-log scale over a limited range, 
can lead to specious power law fits. Hence, and on the ba- 
sis of specific heat measurements alone (i.e. without con- 
sideration of the 1 /Ti /iSR spin-lattice relaxation rate), it 
would therefore appear that the suggestion of unconven- 
tional excitations in Gd2Sn207 should be challenged by 
the principle of "Ockham's razor" . We are therefore led 
to suggest that the description of the specific heat in 
terms of an anomalous power law, ^ T^, in a reduced 
and intermediate temperature range T G [350, 800] mK 
does not provide a convincing indicator for anomalous 
excitations out of the groimd state of Gd2Sn207. Unlike 
the suggestion made in Ref. [53] on the basis of the tem- 
perature independence of the 1 /Ti muon spin relaxation 
rate below Tc ~ 1 K, we have found a fully gapped spin 
wave spectrum with no density of states at zero energy. 
Hence, at this time, the microscopic origin of the temper- 
ature independence of l/Ti found below in Gd2Sn207 
remains to be understood. 



C. Ground state properties 

In the previous section, we have identified a rela- 
tionship between the second J2 and third J31 NN ex- 
change constants which best reproduce the low temper- 
ature thermodynamic behavior in Gd2Sn207. We now 
investigate the role of quantum fluctuations. Fig. [12] dis- 
plays both the minimum spin wave energy gap A and the 
reduction in the staggered moment AS/S along the line of 
best fit J31 = 0.760 J2 - 0.014| Ji|. This result details the 
complicated relationship between the value of the gap, 
and the stability of the ground state, i.e. a decrease in 
the global spin wave energy gap (which may only occur 
at a single fc-point) docs not immediately trigger an in- 
crease of moderate quantum fluctuations. Indeed, the 
opposite behavior is seen in Fig. [T^] In an exchange cou- 
pled Heisenberg antiferromagnet on a non-Bravais lattice, 
the specifics of all relative energy scales come into play, 
and one must not neglect the effects of weakly dispers- 
ing optical modes. All the results presented here can be 
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FIG. 12: (Color online) The global spin wave energy gap (left 
scale) and the reduction in the staggered moment (right scale) 
along the line in the J2 — J31 plane most closely corresponding 
to the parameter set of Gd2Sn207. 



compiled into an effective schematic phase diagram for 
Gd2Sn207. Fig. [13] depicts the separatrix in the J2 — J31 
plane (solid line) that delineates the limit of stability of 
the fc = PC state against a soft mode characterized by 
wavevector 27r/a(3/4, 3/4, 0), i.e. the X— point. Possi- 
bly relevant values of J2 and J31 for which reaches 
it minimum value are shown as the parametric dashed 
line J31 = 0.760 J2 - 0.014|Ji|. We note that Fig. [H 
shows a zero temperature phase diagram that delineates 
the limit of stability of the A; = Palmer-Chalker ground 
state against an instability at 27r/a(3/4, 3/4, 0). A simi- 
lar phase diagram presented in Ref. |23| (with the signs 
corresponding to FM and AF interactions reversed) gives 
the ordering wave vector of the long range magnetic or- 
dered state that first develops as the system is cooled 
down from the paramagnetic phase. 



V. DISCUSSION 

We have considered a Heisenberg model that includes 
isotropic exchange interactions up to third nearest- 
neighbors, single- ion anisotropy and long-range magnetic 
dipole-dipole interactions to describe the long-range or- 
dered state of the Gd2Sn207 pyrochlore antiferromagnet. 
The ground state of this system, as found by neutron 
scattering experiments, corresponds to the (classical, 
PC) ground state described by Palmer and Chalker— for 
the classical Heisenberg pyrochlore with nearest-neighbor 
antiferromagnetic exchange and long-range dipolar cou- 
plingsi^i^i^. We used a long wavelength (1/S spin 
wave) expansion to describe the low-energy excitations 
about the PC ground state and to calculate the low- 
temperature behavior of the specific heat, C^,, and order 
parameter, m, for this material. 

By fitting the available specific heat data in the low- 
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FIG. 13: (Color online) The J2 — J31 ground state phase 
diagram showing regions with order characterized by zero and 
non-zero wavevectors. The solid black line delineates the limit 
of stability of the Palmer-Chalker ground state. The inset 
tetrahedron is tiled with one of the PC states, and the dashed 
line J31 = O.76OJ2 — 0.014| Ji | corresponds to the values of J2 
and J31 which produced the best fits to experimental data. 



temperature range (0.35 K < T < 0.5 K), we were able to 
procure an estimate of the exchange interactions beyond 
nearest neighbors. We obtained evidence that Gd2Sn207 
is in a region of exchange coupling with large stability 
against quantum fluctuations. Our main result (which 
does not rely on excruciatingly fine-tuned exchange con- 
stants beyond nearest-neighbor) is that the experimental 
temperature range 0.35 K < T < 0.5 K corresponds to 
the upper temperature range below which the thermody- 
namic quantities become thermally activated above an 
excitation gap A ~ 1 K^. In other words, the indepen- 
dently experimentally determined microscopic nearest- 
neighbor exchange (on the basis of DC magnetic sus- 
ceptibility), single- ion anisotropy (on the basis of ESR) 
and dipolar coupling strength already predict a temper- 
ature dependence for Cy in Gd2Sn207 that is in rough 
agreement with the experiment without significant ad- 
justment. 

The excitation gap takes its origin from the com- 
bination of single-ion anisotropy and magnetic dipolar 
anisotropy. From our fits of the experimental specific 
heat, we tentatively conclude that the real gap is ac- 
tually even larger than the one we have determined. 
Specifically, considering the lower temperature range in 
Fig. [9] and Fig. [TT] (and the inset of Fig. [Tl] for to), it 
appears that the specific heat is dropping faster in the 
lower temperature range than the calculations predict. 
We speculate that this may indicate that the sub-leading 
anisotropy terms neglected in Ticf in Eq. (5) (and which 
correspond to crystal field terms i?;,,„ with I = 4, 6) 
would further increase the effective gap. In particular. 
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those corrections would resign to further hmit the spin 
fluctuations perpendicular to the local three- fold axis. 
However, at this time, experimental measurements of Cy 
below 0.3 K are required to ascertain quantitatively the 
detail of the microscopic parameters for Gd2Sn207 and to 
determine with better precision the exchange parameters 
Ji, J2 and J31. As in other Gd^+-based insulating mag- 
netic materials^i^, it is also possible that anisotropic 
exchange interactions ultimately need to be included in 
a complete description of Gd2Sn2 07. 

It therefore appears that a rejoinder to the ques- 
tion posed in the Introduction is that the low tem- 
perature specific heat observed in gadolinium stanate 
(Gd2Sn207) may possibly be well described by the con- 
ventional gapped spin wave excitations of Ref . [s^ . We 
believe that cither a confirmation or rebuttal of our 
suggestion of gapped excitations in Gd2Sn207 via spe- 
cific heat (Cy) measurements down to '--^ 100 mK would 
tremendously help focus the discussion about the perva- 
sive low energy excitations in insulating magnetic rare- 
earth pyrochlore oxides. However, a possible confirma- 
tion of such gapped excitations in Gd2Sn207 via specific 
heat measurements would ultimately have to be rational- 



ized within the context of the perplexing and persistent 
temperature-independent spin dynamics found in muon 
spin relaxation studies on this, and other geometrically 
frustrated pyrochlores. 
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